Edge reconstruction in the fractional quantum Hall regime 
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The interplay of electron-electron interaction and confining potential can lead to the reconstruction 
of fractional quantum Hall edges. We have performed exact diagonalization studies on microscopic 
models of fractional quantum Hall liquids, in finite size systems with disk geometry, and found 
numerical evidence of edge reconstruction under rather general conditions. In the present work we 
have taken into account effects like layer thickness and Landau level mixing, which are found to 
be of quantitative importance in edge physics. Due to edge reconstruction, additional nonchiral 
edge modes arise for both incompressible and compressible states. These additional modes couple 
to electromagnetic fields and thus can be detected in microwave conductivity measurements. They 
are also expected to affect the exponent of electron Green's function, which has been measured in 
tunneling experiments. We have studied in this work the electric dipole spectral function that is 
directly related to the microwave conductivity measurement. Our results are consistent with the 
enhanced microwave conductivity observed in experiments performed on samples with an array of 
antidots at low temperatures, and its suppression at higher temperatures. We also discuss the effects 
of the edge reconstruction on the single electron spectral function at the edge. 

I. INTRODUCTION 

The edge of a quantum Hall (QH) system provides a special environment to study electron correlations in one 
dimension. Due to the presence of a strong magnetic field and electron-electron interaction, the bulk of a QH liquid 
is incompressible, while low-lying excitations exist only at the boundary of the liquid. In an experimental sample, 
the physics of edge excitations is strongly affected by the interplay of electron-electron interaction and the confining 
potential due to positive background charge. For example, the edge of an integer QH liquid with a sharp confining 
potential is described by the chiral Fermi-liquid theory and only a single branch of gapless edge excitations exists 
due to the presence of magnetic fieldi When the confining potential is sufficiently smooth, the edge undergoes a 
reconstruction transition in which a portion of the electron liquid is expelled a few magnetic lengths away from the 
periphery of the main droplet. 2,3 Additional low-lying edge excitations that propagate in both directions arise after 
the edge reconstruction transition^ 

The edge excitations of fractional QH liquids are proposed to be described by the chiral Luttinger liquid (CLL) 
theory^ For principal Landau level (LL) filling fractions v — 1/m, the theoretical picture involves only one chiral 
boson mode. On the other hand our exact diagonalization study£ of a microscopic model of fractional QH liquids, in 
finite size systems with disk geometry, suggests that a fractional QH liquid can undergo edge reconstruction for both 
smooth and sharp confining potentials. As a consequence, additional low-energy edge excitations not described by 
the CLL theory are generated and introduce complications into the edge physics. These excitations are clearly visible 
in the low-energy excitation spectrum of a fractional QH system with reconstructed edge. 5 ' 

Apart from purely theoretical interest and importance, our study of fractional quantum Hall edge reconstruction 
has also been motivated by two types of experimental studies of edge physics, both with puzzling results. First, 
the CLL theory predicts a power-law current-voltage dependence (/ ~ V a ) in the tunneling between a Fermi liquid 
metal and a QH edge. The prediction has been tested by experiments&L&Si using samples made by the cleaved- 
edge overgrowth technique^ For a simple filling fraction like v = 1/3, these experiments^*^ found non-Ohmic I-V 
dependence with the exponent a scattering between 2.5 and 2.8, which is close to but noticeably different from the 
CLL theory prediction of a universal exponent a = 3. Furthermore, no convincing plateau behavior away from 
v = 1/3 is presentjL& as predicted by the theory. In fact, for 1/3 < v < 1, a seems to vary continuously for both 
compressible and incompressible values of v, and no universality in a can be extracted from data available to dateiS 
Experimental findings have prompted a number of theories jiLiSiiiiiiiSiiSiiLISi most of which address the continuous 
dependence of a ~ \/v found in one experiment^ only. 

The second type of experiments measure the microwave conductivity of a 2DEG with an array of antidots (mi- 
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croscopic artificial regimes which electrons are forbidden to enter) The microwave conductivity is enhanced with 
a broad peak centered around v = 1/2, exceeding its dc- value by as much as a factor of 5 for microwave frequencies 
up to 10 GHz. The enhanced conductivity is suppressed by increasing temperature and disappears for temperature 
T > 0.5 K. Since the conductivity enhancement and the associated peak are absent for samples without antidots, 
the effect is believed to be related to the antidot edge excitations coupled to the microwave field. However, from the 
submicron size of the antidots (with depletion regime diameter 150 - 250 nm) one can easily estimate the usual edge 
magnetoplasmon modes to have characteristic frequencies above 100 GHz, much higher than the microwave frequen- 
cies used experimentally. Therefore, additional lower-energy edge excitations are required to interpret the mysterious 
data. 

As we discuss later in the paper, it is likely that both sets of puzzles are related to edge reconstruction and in 
particular the additional edge modes associated with it. We present further numerical evidence in support of the 
generality of edge reconstruction in the fractional quantum Hall regime by taking into account effects like layer 
thickness and Landau level mixing, which are of quantitative importance in edge physics but not considered in our 
earlier work. We estimate the finite size effects of our numerical studies and find that they have very little effect 
on our quantitative results of the position of reconstruction transition point. We have also studied in this work the 
electric dipole spectral function that is directly related to the microwave conductivity measurement and the single 
electron spectral function, which is what the tunneling experiments measure. Semiquantitative comparison between 
our results and experiments will be made. 

The rest of the paper is organized in the following way. We review our microscopic model in Sec.[H] Using heuristic 
electrostatic calculations, we discuss the origin of the reconstruction of quantum Hall edges and estimate the finite-size 
effects on microscopic model calculations. In Sec. 11111 we present the numerical evidence for edge reconstruction for 
boundary conditions describing both sharp and smooth edges. We consider the complications due to LL mixing and 
finite thickness of the quasi-2D electron layer. Sec. IIVI studies the effects of thermal fluctuations at finite temperature 
on edge reconstruction. The relevance of our model and the experimentally observed enhanced microwave conductivity 
of samples with antidots is discussed in Sec. El where we present calculations of electric dipole spectral functions. We 
discuss the effects of the edge reconstruction on the single electron spectral function at the edge of fractional quantum 
Hall liquids in Sec. IVII We summarize our results in Sec. IVIII 



We consider a two-dimensional electron gas (2DEG) with disk geometry, as depicted in Fig. da). To model a 
realistic confining potential, as in a modulation-doped AlGaAs/GaAs heterostructure, we assume the neutralizing 
background charge is distributed uniformly on a parallel disk at a distance d above the 2DEG. The radius a of the 
positive charge background is so determined that the disk encloses exactly N / v magnetic flux quanta for ./V electrons 
of the 2DEG, for any desired filling factor v. The bare Coulomb interaction between the background charge and the 
2DEG gives rise to the confining potential. 

Therefore we consider the following Hamiltonian, which describes electrons confined to the lowest LL, using the 
symmetric gauge: 



where cj„ is the electron creation operator for the lowest LL single electron state with angular momentum m. V} nn is 
the matrix element of Coulomb interaction, 



explicitly given by Girvin and JachSi for symmetric gauged U m is the matrix element of the rotationally invariant 
confining potential due to the positive background charge, 



II. MODEL AND ELECTROSTATIC CONSIDERATIONS 




(1) 




(2) 




(3) 



Here </> m is the lowest LL wave function 



4> m {z) = (27r2 m m!)- 1 / 2 z m e-l z l 2 / 4 , 



(4) 
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where z = x + iy is the complex coordinate in the plane of the 2DEG. 

Before we diagonalize the Hamiltonian for microscopic systems to look for evidence of edge reconstruction, we 
first present a heuristic argument that reveals the electrostatic origin of edge reconstruction. Assuming the electrons 
are confined by a hard-wall boundary condition at r = a, we can view the system as a parallel disk capacitor in the 
electrostatic context, as illustrated in Fig.^b). If there is no edge reconstruction, to a good approximation the electron 
density is uniform from the center of the disk all the way to the edge. Together with the uniform positive background 
charge, the capacitor is uniformly charged in this case. Within the electron gas layer, the electrostatic potential is 
a constant in the bulk, but a gradient (or fringe electric field with in-plane component) develops at the edge. This 
fringe field tends to pull the electrons toward the edge; the distance from the edge over which the fringe field effects 
are significant is roughly d. Therefore, the system can gain electrostatic energy by moving electrons outward near 
the boundary along the radial direction. This is expected to happen when the separation between the two oppositely 
charged disks is large enough. When d is large, the fringe field effects are strong; thus the electrostatic energy gain 
from moving some electron density outward at the edge overcomes the associated loss of exchange-correlation energy. 
We identify this as the driving force of edge reconstruction. 

To obtain an estimate of the energetics of this effect, we have calculated the potential energy change of moving one 
electron from one magnetic length Is = y^hc/eB (the typical length scale associated with edge reconstruction) inside 
the edge to the very edge. The energy gain is 

AE = (2d/l B ) tan-\l B /d) + ln(l + d 2 /l%), (5) 

for a half-infinite system, where AE is in units of ve 2 /cIb- This energy diverges as 21n(d/ls) at large d, as shown 
in Fig. |2I a). Once this energy gain exceeds the loss of exchange-correlation energy (which must saturate in the large 
d limit), edge reconstruction occurs. This calculation demonstrates that edge reconstruction must occur in infinite 
systems for sufficiently large d as long as there are states available for electron rearrangements. In addition, it provides 
an estimate of the energy scale associated with edge reconstruction. It also points to a fundamental difference between 
the integer and fractional bulks, in the limit of infinite Landau level spacing. For the former, all the possible single 
electron states are occupied when an infinitely sharp edge is present; thus no reconstruction is possible no matter how 
large d is. On the other hand, this is not the case for fractional bulk and reconstruction is guaranteed to occur for 
sufficiently large d, despite the presence of a sharp edge boundary. 

In this work we perform finite size numerical studies on systems with sizes ranging from 4 to 9 electrons. It is 
important to ask whether the numerical results presented in this paper reveal the physics in the thermodynamic limit 
or merely finite-size artifacts. Repeating the above calculation for a finite-size disk capacitor can provide a measure of 
finite-size effects. The results for systems corresponding to 4 to 9 electrons are shown in Fig. H£b) for d — 21b (which 
is close to critical d c above which edge reconstruction occurs, say at v = 1/3). First, we find that the electrostatic 
energy difference between the edge and one magnetic length inside the edge is very close to the infinite-size value 
(all within 2%). This suggests the finite-size effect is weak, consistent with our finding that there is essentially no 
dependence of d c on size£ Second, finite-size effects reduce the electrostatic energy gain and thus work against the 
edge reconstruction. Therefore, we believe that the edge reconstruction we find is robust and not due to finite-size 
artifacts. 



III. EXACT DIAGONALIZATION STUDY 

In this section, we present the ground state properties and low energy excitation spectra of the model Hamiltonian, 
obtained by exact diagonalization, and look for evidence of edge reconstruction in them. We discuss the appropriate 
choice of boundary conditions at the edge of the system and show that edge reconstruction is a robust property of the 
system under different boundary conditions. In particular, we consider the effects of Landau level mixing near the 
edge due to the presence of a hard-wall potential. We have found that the LL mixing effects, while not of qualitative 
importance to the edge reconstruction physics, do affect the critical spacing d c quantitatively. Our results are also 
robust in the presence of finite thickness of an electron layer, whose effects turn out to be negligible as long as we 
treat the electron layer as a sheet of charge at its maximum density. 

A. Sharp edge 

In the previous study^ we restricted electrons to N/v orbitals of the lowest LL (i.e., no Landau level mixing), from 
in = to m max — N/v — 1. Such a constraint was introduced to describe, for filling factor z/, the presence of a sharp 
cleaved edgeiS beyond which electrons cannot move. In the exact diagonalization study, we found, for v = 1/3, the 
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following evidence that supports the edge reconstruction scenario. (1) The total angular momentum of the global 
ground state M tot becomes greater than that of the corresponding Laughlin state when the separation d between 
electron and background charge layers exceed critical d c = 1.5 ± 0.1 is ; the increase of M to t is due to the outward 
motion of electrons near the edge triggered by the reconstruction, as discussed in the previous section. The critical 
value d c is essentially the same for N = 4-9 electrons, with very little size dependence. (2) The electron density 
profile for the global ground state shows significant oscillation near the edge for d > d c . (3) Counterpropagating 
low- lying excitations can be identified in the low-energy excitation spectrum for d > d c . We also found that property 
(2) is present for filling factors 1/3 < v < 2/3, suggesting edge reconstruction is generic for fractional bulk fillings. 
Properties (1) and (3), on the other hand, are unique to 1/3: (1) because of the existence of the Laughlin state as 
the reference state and (3) because the single chiral boson mode of the unreconstructed edge does not, in general, 
extend to other filling factors. So, the mere existence of counter-propagating modes can not be used for general filling 
factors to look for edge reconstruction. Nevertheless, the general trend remains that M tot increases with d for both 
incompressible and compressible filling fractions, indicating that fractional QH edges tend to reconstruct for large 
d. Figure shows the dependence of Mtot on d for 6, 9, and 12 electrons in 18 orbitals, or v = 1/3, 1/2, and 2/3, 
respectively. In all three cases, Mtot starts to jump at d ~ 11b- 

These numerical findings are consistent with the heuristic arguments that we gave in the previous section. One key 
quantity in these calculations is the critical d, which has always been found to be of the order of I b (even for the more 
complicated cases we discuss later in this section). We believe that this is not a coincidence. Within the lowest LL, Is 
is a fundamental length scale which characterizes the size of single electron wave functions. It is thus also the range 
of effective attraction between electrons due to exchange-correlation effects and therefore the length scale associated 
with edge reconstruction. On the other hand, the separation between the electron and background charge layers, d, 
is the range of the electrostatic fringe field near the edge. Therefore, the electrostatic energy gain generally becomes 
comparable to the loss of electron exchange-correlation energy when d is of the order of Ib] we thus expect d c ~ lg. 

Additional insight into the physics of edge reconstruction may be gained by studying how the low-energy excitation 
spectrum of the system evolves as d increases. According to the CLL theory, for simple bulk fillings like 1/3, there 
exists a single branch of chiral bosonic edge modes with linear dispersion in the long wave-length limit in the absence 
of edge reconstruction; these modes describe the propagation of the deformation of the periphery of a QH system. This 
branch of chiral bosonic edge modes is clearly visible in Figure 01 where we show the low- lying spectrum for an N = 9 
system at v = 1/3 (or, m max — 26), for both d < d c and d > d c . The chirality of the mode is reflected in the fact that 
low- lying excitations exist only in subspaces M > Mq, where Mq is the quantum number of the ground state (which 
is the same as that of the Laughlin state) for d < d c . What is worth emphasizing here are two features not described 
by the CLL theory. First, the low-lying edge modes, which are clearly separated from the bulk excitations with a gap, 
exist for relatively large SM — M — Mq (corresponding to large momentum or short wave-length in translationally 
invariant systems). From the system size we can estimate the momentum k = AM/2wR ps 0.22^ for which the mode 
remains well-defined. Second, the dispersion relation shows significant deviation from linearity and even becomes 
non-monotonous for larger momenta, with a well-defined local minimum near k w 0.152s 1 . This non-linearity grows 
with increasing d. At d — d c , the mode energy crosses zero at this local minimum. Thus the edge reconstruction 
transition can also be understood as an instability, driven by d, of the (unreconstructed) chiral boson mode at finite 
wave vector (of order l^ 1 )- This clearly indicates that edge reconstruction is driven by short-distance physics at the 
scale of Ib, in agreement with previous analysis. Based on this insight, a unified field theoretical description for the 
edges with and without reconstruction, and the transition between the two is advanced by one of us^ 

As discussed earlier, a cutoff at m max for single electron states is introduced to model a sharp cleaved edge. Such 
a sharp edge is not present for all samples. To describe a smooth edge (where confinement is provided only by the 
positive background), we can move the cut-off to higher angular momentum. In fact, even a sharp cleaved edge is 
more complicated than imposing a simple cut-off at m max . More appropriately, electrons are confined in a disk of 
radius a, with a very high step-potential at the edge r = a. The confinement has two related effects. First, the 
step-potential mixes in high Landau level components for single electron wave functions near the edge and raises the 
energies for single-particle states, even with angular momentum m < m max (this is in addition to the effect of the 
back ground charge). Second, the change of single-particle wave function, caused by the in- mixing of LLs, changes 
essentially all the Coulomb interaction matrix elements in Eq. (|TJ . These effects are to be considered in the following 
subsections. 



B. Smooth edge 

As mentioned in the previous subsection, while a sharp cutoff (at m max ) where the background charge ends describes 
samples with cleaved edges, it is not present for other samples. Not surprisingly, we find that removing this sharp 
cutoff can only further favor edge reconstruction. For the 6-electron system considered before^ we have increased 
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the number of orbitals to 30 (m max = 29), while keeping the background charge unchanged so that the sharp cutoff 
is moved about 21b away from the boundary of the background charge. 23 Figures EJa)-(d) show electron densities 
of the global ground states for d = 0.1, 1.6, 2.0, and 5.0 in the presence of such a smoother edge. Compared with 
Fig. 1(b)- (e) in Ref. 0, we find that the change of the sharp cutoff has essentially no effect on the electron density 
for d < 2.0. In particular, the critical distance d c — 1.6 remains very close to that of the sharp edge case. In these 
cases, we expect that the electronic states are almost entirely determined by the competition of the electron-electron 
interaction and the confining potential arising from the background charge. However, with the cutoff moved farther 
out, the edge piece can shift further away from the disk center for larger d so that the total momentum M to t of the 
lowest ground state increases, e.g. from M tot = 65 to 105 for d = 5.0. Therefore, the inclusion of extra single electron 
orbitals has very little effect on the edge reconstruction transition; the effects are important only for systems well in 
the reconstructed phase with d S> d c . For simplicity, in the following discussion we do not include these additional 
orbitals beyond the boundary of the background charge distribution unless otherwise specified. 



C. Layer thickness 



The electrons confined at the interface of the modulation-doped Al^Gai^As/GaAs heterostructures, as used in 
experiments, are not ideally two-dimensional (2D) although their motions perpendicular to the interface arc essentially 
frozen in their ground state due to sharp interface potential. In principle, one needs to use the self-consistently 
calculated wave function appropriate for the heterostructures to study the softening of electron-electron interaction 
due to finite electron layer thickness. However, Stern and Das Sarma 3 ^ showed that the Fang-Howard variational 
wave function^ 

Z (z) = 2{2b)- 3 / 2 ze- z / 2b , (6) 

is a very good approximation to the numerical self-consistent ground state. The Fang-Howard wave function Zq(z) 
peaks at zq = 2b. The parameter b gives the scale of the layer thickness, which is typically ~ 50 A. The finite layer 
thickness weakens the electron-electron interaction, as well as the background charge confining potential. For instance, 
the effective electron-electron interaction in the quasi-2D system is approximated by 

vn? -n e ' Ha a |Zo(*i)H3)(*2)| 2 (7 , 
V (\n -ri\) = ~JJ dz,d Z2 [r2 + {zi _ z2)2]1/2 , (7) 

where r is the in-plane distance between two electrons. 

For typical experimental parameters, the electron wave functions have a finite thickness 2b ~ 100 A, equal to roughly 
one magnetic length. For an AlGaAs/GaAs heterojunction, we choose the z-direction pointing from the AlGaAs to 
the GaAs with the AlGaAs/GaAs interface placed right at the z — plane. Since the AlGaAs introduces a potential 
barrier, the quasi-2D electron density spreads essentially in the GaAs. The background charge layer, introduced by 
(5-doping, is thus located at the z — —d plane (on the AlGaAs side). Experimentally, d ~ 800 A or above. One 
central question, therefore, is whether the edge reconstruction persists in the presence of the finite layer thickness; 
in particular, whether d c for the edge reconstruction remains smaller than the typical d used in experiments. To 
answer this, we repeat exact diagonalization calculations for several values of b, searching for the critical d above 
which the global ground state has a larger total angular momentum than that of the corresponding Laughlin state. 
In this calculation, we use a sharp cutoff (m max = N/u — 1) in the angular momentum space and do not consider the 
complication of LL mixing. Figure summarizes the results for N = 4-9 electrons for v = 1/3. Again, d c is almost 
size-independent for N = 4-9 electrons, confirming that finite-size effects are weak in this calculation. Overall, d c 
decreases as b increases, and can be roughly fit by 

do = d° c - 2b, (8) 

where d° c = 1.4 ± 0.1. Note d c decreases slightly faster than 2b for small b; thus d° c is smaller than d c — 1.5 ± 0.1, 
which is found for systems with zero thickness. Since 2b is the distance between the AlGaAs/GaAs interface and the 
peak of the Fang-Howard wave function, an alternative interpretation of the results is that d c remains roughly as a 
constant, if we measure d from the background charge layer to the plane of maximum electron density, instead of the 
interface. Thus the overall effect of finite layer thickness favors edge reconstruction, and reduces the critical dopant 
layer distance d c slightly. 
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N 


5 6 


7 


8 9 


d c 


0.4 1.6 


0.6 


0.9 1.0 



TABLE I: d c for edge reconstruction for A = (e 2 / 'els) /hu c = 1.0. We define d c as the distance at which Mtot of the global 
ground state exceeds the corresponding value in the d — > limit. 

D. Hard- wall confinement and Landau level mixing near the edge 

So far, we have been working in the limit where the kinetic energy is completely quenched, so electrons are entirely 
in the lowest LL. In this limit, the Coulomb energy e 2 /els is the only energy scale in the system. However, for 
typical n-type GaAs heterostructures in experiments, the Coulomb energy is comparable to the cyclotron energy huj c 
separating the LLs. For the fractional QH effect, the LL mixing effects have been considered in numerical as well as 
analytical studies. 24 i 25 i 26 i 27 

The LL mixing effects are important for edge physics, since the single particle energy for electrons confined in a 
disk increases monotonically from bulk to edge and eventually crosses higher LL energies due to confinement, as first 
discussed by HalperinSS in the QH context. Therefore, LL mixing, in particular, resulting from the cleaved sample 
edge may as well alter the edge physics. To include this effect, we solved the Schrodinger equation for non-interacting 
electrons confined in a two-dimensional disk (we neglect the finite thickness of the electron layer) by a hard-wall 
boundary condition: 

ip(r) = 0, for r > a. (9) 

The ground state wave function (r) in each angular momentum (m) subspace now becomes a mixture of states 
in all LLs with the same to quantum number. As a result, the energy of the ground states increases from %ui c /2 
(lowest LL value) for to = to approach 3Tuu c /2 (first LL value) for m max = 3N — 1, as depicted in Fig. [7{b) 
(here we do not include contributions from the confining potential of the background charge). Figure E{ a) shows 
the cumulative overlaps, Y27=o |(V'm( r )l < / ) m( r ))| > °f the corresponding ground state wave functions [V'mM] f° r eacn 
angular momentum to, with the LL wave function [<^n( r )] f° r the lowest five LLs (i = 0-4). The total overlap of the 
five lowest LLs is more than 99% for each m. 

Following the usual procedure of projecting the Hamiltonian of the system onto the ground state manifold, we 
obtain an effective Hamiltonian in the same form as Eq. ^ with cj^ creating an electron in the LL-mixed ground 
state with angular momentum to. The single particle wave functions cf) m in Eq. (J2J and © are replaced by ip!^. We 
absorb the m-dependent single-particle energy e m into the confining energy term U m . We restrict the Hilbert space to 
< to < 3N for v — 1/3. Note that we have, in addition to e 2 /ds, a second energy scale Tiuj c . The LL mixing effect 
is thus characterized by the dimensionless parameter A = (e 2 /eZs)/?iw c . We are, in particular, interested in A ~ 1, 
which is close to the real experimental conditions. The LL mixing raises e m by roughly hu c for edge states, thereby 
making electrons occupying these states energetically unfavorable. However, the squeezing of the wave functions near 
the edge due to the Landau level mixing effect lowers the confining potential from the background charge and also 
reduces the range of effective attraction due to exchange-correlation effects. These effects favor edge reconstruction. 

Figure Q3] shows the low-energy spectra for N = 9 electrons in 27 orbitals with hard- wall boundary conditions for 
various d and A = 2.0. After edge reconstruction, the total ground state momentum becomes M to t = 116, increasing 
from M to t — 108 of the corresponding Laughlin state. Counterpropagating low-energy modes can be observed near 
Mtot — 116 for d — 1Mb- Qualitatively, these results agree very well with the scenario using sharp cutoff in angular 
momentum space (Fig. [1J. However, we also observe fluctuations in the critical d c for systems with different sizes. 
For A = 1.0, we list d c in Tabled which varies from 0.4 Ib to 1.6 Ib- Although the data seems to stabilize at d c w 1.0, 
we cannot draw definitive quantitative conclusions without data from larger systems in contrast to the case of no LL 
mixing. We believe the increased finite size effect here is due to the fact that LL mixing affects states in a relatively 
wide region near the edge; this effect was not present in our earlier study. We also point out that here we define 
d c as the distance at which M to t exceeds the corresponding value in the d — ► limit and the results presented are 
based on this working definition. With strong LL mixing effects, the ground state for d < d c may not have an M to t 
consistent with the value of the corresponding Laughlin state. We believe that such a different M tot in the limit of 
d — > suggests the Laughlin state may not be a good approximation to the ground state. In other words, we probably 
cannot identify the finite systems as having filling fraction 1/3 unambiguously. 

We would like to emphasize that all of our numerical results, as well as the heuristic arguments based on electrostatic 
considerations, suggest that d c ~ Ib- In real samples, on the other hand, d > 101b- We thus believe that it is safe to 
conclude that the edges of all samples studied so far are reconstructed when the bulk filling is fractional, regardless 
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of whether the edges are cleaved or not. 

IV. EFFECTS OF FINITE TEMPERATURE ON EDGE RECONSTRUCTION 

Thus far our studies have been focusing on ground state (or zero temperature) properties of the system, especially 
those related to edge reconstruction. It is of interest to investigate how finite temperature, and the thermal fluctu- 
ations associated with it, affect these properties. In our previous study^ we used a finite-temperature Hartree-Fock 
approximation to show that the reconstruction of a v = 1 QH edge is suppressed above a certain temperature T* . This 
is quite reasonable since thermal fluctuations tend to suppress edge density oscillations associated with reconstruction. 
In that study the temperature T* m 0.05e 2 /cIb (or T* s=s 6 K) for typical experimental parameters. It is expected, 
however, that the Hartree-Fock calculation tends to overestimate the temperature scale, due to effects of finite size, 
layer thickness, LL mixing, disorder, etc. Here we study the effects of finite T for fractional bulk filling v = 1/3. In 
this case a Hartree-Fock calculation is no longer possible; we thus use exact diagonalization to generate all the excited 
states with significant thermal weight and perform thermal averaging. Figure shows the evolution of the density 
profile for 7 electrons at v — 1/3 at several temperatures. We find that, similar to what we found earlier for v = 1, the 
edge density oscillation associated with edge reconstruction is washed out above a certain characteristic temperature, 
T* w 0.05e 2 /eZs, at which the density profile no longer has strong oscillations and becomes quite similar to that of 
the corresponding Laughlin state. This is about the same temperature scale (but for much smaller systems) that we 
got for the integer edge from Hartree-Fock; we expect the temperature scale will become much lower for a comparable 
size system. 

V. MICROWAVE CONDUCTIVITY AND ITS TEMPERATURE DEPENDENCE 

The microwave conductivity related to microwave absorption of our microscopic system can be calculated through 
the electric dipole spectral function as follows: 

a(u) « l^/le • ^}| 2 <5 (Ef - Ek - fiu) f(Ek) (E f - E. t ) , (10) 
»-►/ 

where f(E) — e~ E ' ksT /Z is the thermal weight of a state with energy E. Eij and ipij are the energies and wave 
functions for initial and final states, respectively, and the summations are over all eigenstates. e is the unit vector of 
an electric field and r the position vector for electrons. 

The microwave conductivity, as determined from Eq. (|10|) for a system of 6 electrons for v = 1/3 at different 
temperatures, both before id = 1.01b) and after edge reconstruction (d — 2.01b), are shown in Figs. HOf a) and (b). 
The data is coarse-grain averaged in frequency space for clarity of viewing. For d = 1Mb, the spectral function 
shows a dominant peak at u> » 0.04, in units of e 2 /hels- This peak corresponds to the single chiral edge mode for the 
principal filling factor v — 1/3, marked by a solid arrow in the corresponding excitation spectrum, Fig JlOf c). At higher 
temperature, the peak becomes less prominent, due to the reduced statistical weight of the ground state and other 
low-lying states; however the position of the peak does not shift with changing T, and the peak can still be identified 
as the dominant one in the whole spectrum for temperature up to T — 0.05, in units of e 2 /els- For d — 2.01b, where 
edge reconstruction has occurred, the frequency dependence of the conductivity becomes qualitatively different. Two 
distinct peaks can be resolved for lu < 0.05 at T = 0.01, contributed by three dominant modes, with two in the lower 
peak (marked by a dotted arrow and a dashed arrow). The additional modes are due to edge reconstruction, which 
creates two counter-propagating edge modes. What is more interesting is the manner in which the conductivity evolves 
as T increases; here we find, at higher temperatures, the peak of the spectral function shifts to higher frequencies, 
and the low frequency response due to the additional modes get suppressed. The low-temperature peaks can barely 
be resolved at T = 0.05. Such behavior is in good qualitative agreement with the microwave experiments in samples 
with an array of antidots (reported in Ref. 1 1 9^) . These authors find enhanced low frequency (much lower than the 
edge magneto plasmon frequency) conductivity at low T, while such enhancement gets suppressed at higher T . One 
should note, however, that the system size of our study (as, say, parameterized by the circumference of the edge) is 
much smaller than those of the real samples; thus, not surprisingly, the energy and temperature scales obtained in 
our work is considerably larger than those of the experimental data; this is purely a finite size effect. On the other 
hand the temperature scale obtained here is consistent with what we obtained in the previous section through the 
temperature dependence of edge density oscillation. 

The microwave experiment on samples with antidots 19 > suggests that the observed enhanced conductivity is a generic 
feature for all fractional filling factors and peaks around v = 1/2. Unfortunately, finite-size effects do not make it 



8 



possible for us to obtain conclusive results for arbitrary v, except for simple cases such as v = 1/3. This is because the 
hierarchy FQH states have more complicated edge structure, which necessarily leads to stronger finite-size effects; the 
situation is even worse for compressible bulk like that that of v — 1/2, because the gapless bulk and edge excitations 
are inevitably mixed together. Here we discuss another relatively simpler case, v = 2/3, where the QH liquids can 
be understood as a v — 1/3 hole condensate embedded in a v = 1 electron condensatei^LiSSSi Therefore, two 
counterpropagating edge modes exist even without the edge reconstruction. We find, not surprisingly, these two 
modes (AM = M to t — Mo — ±1) in °ur numerical calculation for d — 0.1, as shown in Figs. lllf a') and (c). One mode 
(AM = 1, marked by a solid arrow) has a squared dipole matrix element almost two orders of magnitude larger than 
the other (AM = — 1, marked by a dotted arrow), so the low-temperature spectral function is dominated by a single 
mode (the edge magnetoplasmon mode). On the other hand, for d = 2.0, we find, more edge modes which include two 
AM = — 1 low-energy excitations as shown in Fig. II If . These new modes can roughly be regarded as the results 
of the edge reconstruction of the inner v = 1/3 hole condensate. There is no evidence that the outer v = 1 edge can 
be reconstructed since we model samples with sharp cleaved edges. The dipole spectral function, however, is now 
dominated by one of the AM = — 1 modes [marked by a solid arrow in Fig. Illf c)]. with a squared spectral weight 
more than two orders of magnitude greater than those of the AM = 1 modes (the lowest one marked by a dotted 
arrow). Thus, quite similar to the v = 1/3 case, we find additional modes due to edge reconstruction at larger d. We 
also find the transition processes that dominate the dipole spectral function and microwave conductivity change due 
to edge reconstruction. Also similar to the v = 1/3 case, the dipole spectral function becomes dominated by bulk 
excitations at high temperatures and the low-frequency spectral weight gets suppressed. Thus the calculation of the 
dipole spectral function at v = 2/3 also finds additional low-energy modes generically arise from the reconstruction 
of fractional QH edges, and they make important contribution to the low-frequency dipole spectral function. Such 
contribution, however, get suppressed at high temperatures. Thus this behavior is not specific to principal bulk filling, 
in agreement with the enhanced microwave conductivity in antidot samples in the entire fractional filling range, and 
its suppression at higher temperature. 

VI. SINGLE ELECTRON SPECTRAL FUNCTION AND EDGE TUNNELING 

Numerical calculations of single electron spectral function at the edge in finite size systems have been performed 
by Palacios and MacDonald^ They considered a QH droplet with Coulomb interaction, but without a physically 
realistic confining potential. In particular, they calculated the squared matrix elements between the ground state of an 
TV-electron system to the low- lying states of the corresponding (N + l)-electron system at v — 1/3. These numerical 
results can be compared to those obtained by the CLL theory, which predicts that the low-lying energy spectrum of 
a QH system at principal filling fraction, such as v — 1/3, can be described by a branch of single-boson edge states 
with angular momentum / (I = 1,2,3, ...) and energy ei. In the CLL language, we can label each low-energy state by 
a set of occupation numbers {n/}, whose total angular momentum and energy are M = Mo + AM = Mo + J2i 
and E — Eq + AE = Eq + X); n ' e 'i respectively, where, Mo and Eq are total angular momentum and energy of 
the corresponding ground state. Palacios and MacDonald 35 found excellent agreement between the squared matrix 
elements, T({ni}) — | (?/>{ ni } (TV + l)\cl N+AM \ifJo(N))\ 2 , calculated numerically in the microscopic model and those 
calculated based on the CLL theory. Note that Mq(N + 1) — Mq(N) = 3N is the difference in total angular momenta 
between the N- and (N + l)-electron ground states. Such a comparison is made possible by the unambiguous 
identification of the low energy spectrum in terms of {ni} (based on AM and AE), as well as in the agreement of the 
corresponding values of T({ni}), at least for M < 4. 

In this section, we study the single electron spectral function in the presence of a physically realistic edge confining 
potential, generated by a layer of background charge, distributed uniformly on a disk at distance d above the 2DEG. 
As discussed earlier, the single electron spectral function is directly measured in edge tunneling experiments. We are 
interested in finding out how the electron spectral function is affected by the presence of the edge confining potential, 
which mixes the eigenstates in the system without the confinement. In particular, we are interested in whether and 
how the electron spectral function is affected by edge reconstruction. 

We calculated the tunneling spectral weights for v — 1/3 by adding one more electron into a system of 6 electrons. 
The realistic edge confining potential is generated by the appropriate background charge, which neutralizes the 
resulting system. For d — 1.0 (before reconstruction), we plot the spectrum and electron spectral weights in Fig.^Ja) 
and (b), respectively. As in the absence of the confining potential, we can identify the lowest edge excitations for 
each M = Mo + AM subspace in the low-energy excitation spectrum as the single-boson edge state with ni = <5/,am 
(with excitation energy ei). Thus, the family of edge states can be unambiguously identified, since their AM and AE 
must be simultaneously written as linear combinations of I and e; , respectively. We list the tunneling spectral weights 
for AM < 4 in Table [HJ which are highlighted by solid lines in Fig. ^| These matrix elements are rather close to 
those obtained in the absence of the confining potential^ and consistent with the predictions of the CLL theory (for 
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{1100} 
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3.644 


4.5 




{0010} 


0.0461 
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0.0306 
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{4000} 


0.1252 


2.083 


0.0972 


3.413 


3.375 




{2100} 


0.1038 


5.182 


0.0797 


5.113 


6.75 




{1010} 


0.0756 


2.529 


0.0535 


4.312 


3 




{0200} 


0.0853 


0.587 


0.0621 


0.955 


1.125 




{0001} 


0.0465 


0.402 


0.0266 


0.839 


0.75 



TABLE II: Tunneling spectral weights T({n;}) = \(ip{ n y(N + 1) |c3 JV+AM IV'o(Af)) | 2 , for microscopic model at v = 1/3, before 
(d — 1.0/b) and after (d = 1.61b) edge reconstruction transition (normalized to the ground state-to-ground state matrix element, 
To), and for CLL theory. The microscopic system contains N = 6 electrons with corresponding background charge before an 
additional electron tunnels into the electron layer. AA1 = M — Mo and AE = E — Eo (in units of e 2 /eIb) are total angular 
momentum and energy measured from the corresponding ground state values (Mo and -Bo) for the resulting system, {n;} is a 
set of occupation numbers of the bosonic edge- wave with angular momentum I and energy e; . 

infinite system). 

For reconstructed edges on the other hand, the situation becomes more complicated since there are additional 
nonchiral boson excitations, all of which are coupled in general. However, it has been proposed^ that in the strong 
coupling limit (for the unscreened long-range Coulomb interaction which has logarithmic singularity in the long-wave 
length limit) one mode, which represents the total charge density mode, may dominate and behaves just like the single 
branch below the edge reconstruction transition. For d = 1.6 (after reconstruction), we plot the tunneling spectrum 
and spectral weights in Fig. H"27 c) and (d), respectively. For the lowest two states with AM = +1, we find that the 
squared tunneling matrix element for the lower state is 0.213 and for the second lowest state is 2.845; this is close 
to 3 which is predicted for the {1000} mode in the CLL theory. For this reason, we may identify the two states as 
members of the neutral mode and the charge mode, respectively (also the charge mode is indeed expected to have 
higher energy). In addition, we find that the corresponding squared matrix elements for the ground states in AM = 
2-4 subspaces are 1.681, 1.320, and 0.839, respectively. These numbers are close to 1.5, 1, and 0.75 predicted for the 
edge modes in the corresponding subspaces. So these three states may be identified as {0100}, {0010}, and {0001} 
of the charge mode. If the distinction of charge and neutral modes can be made in this way and the charge mode 
indeed controls the tunneling behavior, we should be able to generate a family of excitations with angular momenta 
and energies which can be calculated according to the four charge excitations identified so far. We should also be 
able to find, near the calculated energies in the corresponding angular momentum subspace, excitations with squared 
matrix elements close to the predictions in the CLL theory. This indeed seems to be the case, as shown in Table ITT1 
The listed states are, again, highlighted in Fig. ^| We emphasize that the squared matrix elements for the rest of 
the low-lying states (AE < 0.05e 2 /ds) with AM ^ are typically very small, with the largest one being roughly 
unity (or the ground state-to-ground state value) . We point out that in order to calculate the matrix elements for up 
to AM = +4, we use a smooth edge with m max = 23. We found that reducing m max (giving a sharper edge) has no 
significant effects on the squared matrix elements in the reduced subspaces. Our results thus suggest that the effect 
of edge reconstruction on the structure of the single electron spectral function is fairly weak; this is consistent with 
the experimental finding that the tunneling exponent is close to the prediction of the CLL theory at v = 1/3, despite 
the fact that the edges are expected to be reconstructed. 
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VII. CONCLUSIONS 



In this paper, we have performed exact diagonalization studies on microscopic models of fractional QH liquids 
in systems with disk geometry, and investigated the interplay between electron-electron interaction and confining 
potential due to background charge near the edge. We have shown that the edges of fractional QH liquids reconstruct 
when the background charge (dopant) layer is separated far enough from the electron layer, and the critical distance 
for this to happen is of order one magnetic length. The edge reconstruction happens because the electrostatic energy 
gained by moving electrons outward near the sample edge increases logarithmically as the separation increases, and 
eventually exceeds the loss of electron exchange-correlation energy. Such behavior is found to have very weak finite-size 
effects in most cases, even for small systems with 6-9 electrons that we used in our studies. Our results suggest that 
edge reconstruction occurs rather generically in high-quality AlGaAs/GaAs samples used in experimental studies, as 
the corresponding distance in these samples are typically of the order of ten magnetic lengths or larger. 

In our studies we have used different types of boundary conditions for electronic wave functions near the edge, 
corresponding to different types of samples (e.g., whether edge is cleaved or not). We have demonstrated that the 
edge reconstruction phenomenon is not sensitive to the choices of specific boundary conditions qualitatively, be it 
hard-wall confinement or one that leads to a smooth confining potential. While different boundary conditions lead to 
quantitatively different critical spacing between dopant and electron layers, our conclusion that real samples are all 
in the reconstructed regime is robust. 

With reconstructed edges, fractional QH liquids can have additional edge modes that propagate along both direc- 
tions. In general, we find these modes tend to have much lower energy scales than the edge modes in the absence of 
edge reconstruction. Therefore, they can have very important effects on the low-energy behavior of edge transport 
and tunneling experiments. We have performed calculations on the electric dipole spectral function as well as single 
electron spectral function, for systems with and without edge reconstruction. We find that edge reconstruction affects 
the dipole spectral function rather strongly, and its frequency as well as temperature dependences compare favorably 
with microwave conductivity measurements performed in samples with an array of antidots (and their associated 
edges) . On the other hand we find the electron spectral function at v = 1/3 is not modified strongly by edge recon- 
struction; this is consistent with tunneling experiments which find the tunneling exponent at v = 1/3 quantitatively 
close to the prediction of the chiral Luttinger liquid theory despite the presence of edge reconstruction. 
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FIG. 1: (a) Sketch of the system with rotational symmetry considered here, which is made of an electron layer and a uniformly 
distributed, neutralizing background charge layer separated by a distance d from each other. Electrons are confined by a 
hardwall boundary condition, so they cannot move beyond the edge of the background charge, (b) The side view of the system. 
If electrons are uniformly distributed, the electrostatic potential is a constant in the bulk of the electron layer, but a gradient 
(or fringe electric field with in-plane component) develops at the edge, which tends to pull the electrons toward the edge. 
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FIG. 2: Estimate of electrostatic energy associated with edge reconstruction and finite size effects. Consider a parallel 
disk capacitor with uniformly distributed positive and negative charges. Left panel shows the potential energy gain [AE = 
(2d/is) tan _1 (is/d) + ln(l + d 2 /l%)] in moving an electron from one magnetic length (Ib) inside the edge to the edge, for a 
half-infinite capacitor system (the infinite-size limit of the disk system) . d is the distance between the two charge layers. Right 
panel shows AE in finite-size disk systems (4-9 electrons), with d = 21b, at v — 1/3. The dashed line is the infinite-size limit 
for AE at d = 2l B . 
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FIG. 4: Low energy spectra for N — 9 electrons in 27 orbitals for (a) d = 0.5, (b) d = 1.0, (c) d = 1.5, and (d) d = 2.0, in 
units of Ib- Excitation energies (A_E) are measured, in units of e /e/s, from the ground state in the Mtot = 108 (that of the 
corresponding Laughlin state) subspace. The edge reconstruction transition occurs around d = 1.5, as AE becomes negative 
for M tot = 115. 




FIG. 5: The electron density p(r) of the global ground state for 6 electrons in 30 orbitals compared with that of the Laughlin 
state (dotted lines) for (a) d = 0.1, (b) d = 1.6, (c) d = 2.0, and (d) d = 5.0, in units of Ib- The radius of the disk with uniform 
background charge corresponds to v — 1/3, so electrons are allowed to move ~ 21b beyond the background charge. 
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FIG. 6: Critical d c (in units of Ib) for edge reconstruction, at which the ground state momentum Mtot becomes greater than 
that for the corresponding Laughlin state, for N = 4-9 electrons at v = 1/3 with finite layer thickness. d c is measured as 
the distance from the positive background charge layer to the (GaAs/AlGaAs) interface where potential is discontinuous. The 
finite thickness of the 2D electron gas in the perpendicular direction is described by the Fang-Howard variational wave function 
Z (z) = 2{2b)-' i/2 ze- z/2b (see inset). d c can be roughly fit to d c = 1.4 — 2b, where 1.4 ± 0.1 can be regarded as the critical d 
for zero layer thickness, and 2b the distance from the peak of the variational wave function to the interface (z = 0). 
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FIG. 7: (a) Cumulative overlaps, y^"_ n | (Ol^mM) | , of the ground state wave function [ip^ (r)] for each angular momen- 
tum m, in the presence of a hard- wall boundary condition, with LL wave function, 4> % m {r), for the lowest five LLs (i — 0-4). 
The sum of the overlaps is more than 99% for each m. (b) The single-particle energy of the ground state e m increases from 
TiWc/2 in the bulk to 3hu> c /2 at the edge. 
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FIG. 8: The low energy spectra for N — 9 electrons in 27 orbitals with hard- wall boundary conditions for (a) d — 0.51b, 
with no edge reconstruction, (b) d = 1.01b, close to but before edge reconstruction, (c) d = 1Mb, and (d) d = 2.01b, both 
after edge reconstruction. Excitation energies (AE) are measured, in units of e 2 /eZs, We choose the dimensionless parameter 
A = (e 2 /elB)/Tujj c — 2.0 here. After edge reconstruction, the total ground state momentum becomes M to t = 116, increasing 
from Mtot = 108 of the corresponding Laughlin state. 
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FIG. 9: Finite-temperature (in unit of e 2 /ds) electron density profiles p(r) for N — 7 electrons at v = 1/3. The distance 
between the two charge layers is fixed at d = 31b- At and above T = 0.05e /(1b, p{r) become similar to that of the 
corresponding Laughlin state, fluctuating slightly around 1/3. This is very different from p(r) at T = 0.005e 2 /d_g, where 
strong density oscillation due to edge reconstruction can be seen. 
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FIG. 10: Microwave conductivity a(ui) calculated by electric dipole spectral function for N = 6 electrons at v — 1/3 for 
(a) d = 1Mb (with no edge reconstruction) and (b) d — 2Mb (after edge reconstruction) at various temperatures. The 
corresponding low-energy excitation spectra are plotted for (c) d = 1Mb and (d) d = 2Mb- Before edge reconstruction, the 
absorption is dominated by long wave-length chiral edge mode [solid arrows in (a) and (c)] even above T = 0.05e 2 /el. After 
reconstruction, extra modes [indicated by three arrows in (b) and (d)] due to reconstructed edge can be identified as contributing 
to the peaks in a{u) for u> < 0.05, in units of e 2 /el. At T > 0.03e 2 /el, these modes become less significant, compared to bulk 
excitation. 
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FIG. 11: Microwave conductivity (j(uj) calculated by electric dipole spectral function for N = 12 electrons at v — 2/3 for 
(a) d = 0.11b and (b) d = 2Mb at various temperatures. The corresponding low-energy excitation spectra are plotted for (c) 
d = 0.11b and (d) d = 2.01b- For d — 0.11b, f(w) is dominated by the lowest AM = +1 edge mode [solid arrows in (a) and (c)], 
while the lowest AM = — 1 edge mode [dotted arrows in (a) and (c)] is significantly weaker. For d — 2.01b, ct{uj) is dominated 
by the second lowest A M = — 1 edge mode [solid arrows in (b) and (d)]. Both the lowest AM = ±1 modes are much weaker, 
as indicated by dashed and dotted arrows in (b) and (d). 
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FIG. 12: (a) Spectrum and (b) normalized spectral weights, T({n;}) /To, for tunneling one additional electron into a 6-electron 
system at v — 1/3 for d — 1.01b- (c) Spectrum and (d) normalized spectral weights, T({n;})/To, for tunneling one additional 
electron into a 6-electron system at v = 1/3 for d = 1Mb- In all figures, solid lines repreent states that have significant 
contribution to the electron spectral function, whose corresponding matrix elements are listed in Table HH 



